use "/winemag-final.dta"

//regressions for POINTSGROUP//
ologit pointsgroup wc sixltr article ipron prep auxverb verb adj compare posemo tentat certain reward discrep social affiliation family leisure percept see hear bio body health ingest, ro
margins, dydx(var) predict(pr outcome(1)) predict(pr outcome(2)) predict(pr outcome(3)) predict(pr outcome(4))
marginsplot, recast(scatter)

//marginal effects//
//closed class terms//
ologit pointsgroup wc sixltr article ipron prep auxverb verb adj compare posemo tentat certain reward discrep social affiliation family leisure percept see hear bio body health ingest, ro

//closed class words//

margins, dydx(wc)
marginsplot, recast(scatter) xlabel(1 "Value" 2 "Popular" 3 "Premium" 4 "Super Premium" 5 "Ultra Premium" 6 "Luxury" 7 "Super Luxury" 8 "Icon")

margins, dydx(sixltr)
marginsplot, recast(scatter) xlabel(1 "Value" 2 "Popular" 3 "Premium" 4 "Super Premium" 5 "Ultra Premium" 6 "Luxury" 7 "Super Luxury" 8 "Icon")

margins, dydx(ipron)
marginsplot, recast(scatter) xlabel(1 "Value" 2 "Popular" 3 "Premium" 4 "Super Premium" 5 "Ultra Premium" 6 "Luxury" 7 "Super Luxury" 8 "Icon")

margins, dydx(auxverb)
marginsplot, recast(scatter) xlabel(1 "Value" 2 "Popular" 3 "Premium" 4 "Super Premium" 5 "Ultra Premium" 6 "Luxury" 7 "Super Luxury" 8 "Icon")

margins, dydx(article)
marginsplot, recast(scatter) xlabel(1 "Value" 2 "Popular" 3 "Premium" 4 "Super Premium" 5 "Ultra Premium" 6 "Luxury" 7 "Super Luxury" 8 "Icon")

margins, dydx(verb)
marginsplot, recast(scatter) xlabel(1 "Value" 2 "Popular" 3 "Premium" 4 "Super Premium" 5 "Ultra Premium" 6 "Luxury" 7 "Super Luxury" 8 "Icon")

margins, dydx(adj)
marginsplot, recast(scatter) xlabel(1 "Value" 2 "Popular" 3 "Premium" 4 "Super Premium" 5 "Ultra Premium" 6 "Luxury" 7 "Super Luxury" 8 "Icon")

margins, dydx(prep)
marginsplot, recast(scatter) xlabel(1 "Value" 2 "Popular" 3 "Premium" 4 "Super Premium" 5 "Ultra Premium" 6 "Luxury" 7 "Super Luxury" 8 "Icon")

graph combine wc-pointsgroup.gph sixltr-pointsgroup.gph ipron-pointsgroup.gph auxverb-pointsgroup.gph articles-pointsgroup.gph verb-pointsgroup.gph adj-pointsgroup.gph prep-pointsgroup.gph

//persuasion//

margins, dydx(compare)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(posemo)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(tentat)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(certain)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(cogproc)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(reward)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(discrep)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

graph combine compare-pointsgroup.gph posemo-pointsgroup.gph reward-pointsgroup.gph tentat-pointsgroup.gph certain-pointsgroup.gph discrep-pointsgroup.gph

//social experience//

margins, dydx(social)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(affiliation)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(family)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(leisure)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

graph combine social-pointsgroup.gph affiliation-pointsgroup.gph family-pointsgroup.gph leisure-pointsgroup.gph

//somatic experience//   

margins, dydx(percept)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(see)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(hear)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(bio)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(body)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(health)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")

margins, dydx(ingest)
marginsplot, recast(scatter) xlabel(1 "Classic" 2 "Outstanding" 3 "Very good" 4 "Good")


graph combine percept-pointsgroup.gph see-pointsgroup.gph hear-pointsgroup.gph bio-pointsgroup.gph body-pointsgroup.gph health-pointsgroup.gph ingest-pointsgroup.gph


value ($4–$10), popular ($10–$15), premium ($15–20), super premium ($20–$30), ultra premium ($30–$50), luxury ($50–$100), super luxury ($100–$200) and icon ($200 or more).

label define pricecode 0 "Value ($4–$10)" 1 "Popular ($10–$15)" 2 "Premium ($15–20)" 3 "Super premium ($20–$30)" 4 "Ultra premium ($30–$50)" 5 "Luxury ($50–$100)" 6 "Super luxury ($100–$200)" 7 "Icon ($200 or more)"
